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Summary. Wc study supercritical branching Brownian motion on the real 
line starting at the origin and with constant drift c. At the point x > 0, we 
add an absorbing barrier, i.e. individuals touching the barrier are instantly 
killed without producing offspring. It is known that there is a critical drift 
Co, such that this process becomes extinct almost surely if and only if c > cq. 
In this case, if denotes the number of individuals absorbed at the barrier, 
we give an asymptotic for P{Zx = n) as n goes to infinity. If c = cq and the 
reproduction is deterministic, this improves upon results of L. Addario-Berry 
and N. Broutin [1] and E. Ai'dckon [2] on a conjecture by David Aldous about 
the total progeny of a branching random walk. The main technique used in 
the proofs is analysis of the generating function of near its singular point 
1, based on classical results on some complex differential equations. 

Keywords. Branching Brownian motion, Galton-Watson process, Briot- 
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generating functions. 
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1 Introduction 

We define branching Brownian motion as follows. Starting with an initial individual sitting 
at the origin of the real line, this individual moves according to a 1-dimensional Brownian 
motion with drift c until an independent exponentially distributed time with rate 1. At that 
moment it dies and produces L (identical) offspring, L being a random variable taking values 
in the non-negative integers with P{L = 1) = 0. Starting from the position at which its 
parent has died, each child repeats this process, all independently of one another and of their 
parent. For a rigorous definition, see for example [10]. 

We assume that m = E[L\ — 1 G (0, oo), which means that the process is supercritical. At 
position X > 0, we add an absorbing barrier, i.e. individuals hitting the barrier are instantly 
killed without producing offspring. Kesten proved [19] that this process becomes extinct 
almost surely if and only if the drift c > cq = \/2m (he actually needed E[L'^] < oo for the 
"only if" part, but we are going to prove that the statement holds in general). A conjecture 
of David Aldous [3], originally stated for the branching random walk, says that the number 
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A'^^^ of individuals that have lived during the lifetime of the process satisfies -E'[A'a;] < oo and 
E[Nx log^ Nx] = oo in the critical speed area (c = cq), and P{Nx > n) ~ Kn~'^ in the 
subcritical speed area (c > cq), with some > 0,7 > 1. For the branching random walk, 
the conjecture of the critical case was proven by Addario-Berry and Broutin [1] for general 
reproduction laws satisfying a mild integrability assumption. Ai'dekon [2] refined the results 
for constant L by showing that there are positive constants p,Ci,C2, such that for every 
X > 0, we have 

— < P{Nx > n) < — for large n. 

n[iogny nyiogny 

Assuming L constant has the advantage that Nx is directly related to the number Zx of 
individuals absorbed at the barrier by Nx — 1 = {Zx — 1){L/{L — 1)), hence it is possible to 
study Nx through Zx- In this sense, Neveu [23] had already proven the critical case conjecture 
for branching Brownian motion since he showed that the process Z = {Zx)x>o is actually a 
continuous-time Galton- Watson process of finite expectation, but with E[Zx log^ Zx] = 00 
for every x > 0, if c = cq. 

Let N = {1, 2, 3, . . .} and No = {0} U N. Define the infinitesimal transition rates (see [4], 
p. 104, Equation (6) or [14], p. 95) 

qn = liin-PiZx=n), n G No\{l}. 

2:4.0 X 

We propose a refinement of Neveu's result: 

Theorem 1.1. Assume c = cq. Assume that £'[L(log L)^^'^] < 00 for some e > 0. Then we 
have as n ^ 00, 

% ~ n "■'^d P{Zx > n) ~ 'J^^ for each x > 0. 

^-^ 711 log nl^ n(iogn)^ 

k=n 

The heavy tail of Zx suggests that its generating function is amenable to singularity 
analysis in the sense of [12]. This is in fact the case in both the critical and subcritical cases if 
we impose a stronger condition upon the offspring distribution and leads to the next theorem. 

Define /(s) = £'[s^] the generating function of the offspring distribution. Denote by b the 
span of L — 1, i.e. the greatest positive integer, such that L — 1 is concentrated on (5Z. Let 
Ac < Ac be the two roots of the quadratic equation A^ — 2cA + Cq = and denote by d = ^ 
the ratio of the two roots. Note that c = cq if and only if Ac = Ac if and only if d = 1. 

Theorem 1.2. Assume that the law of L admits exponential moments, i.e. that the radius 
of convergence of the power series E[s^] is greater than 1. 

— In the critical speed area (c = cq), as n ^ 00, 

Cq CqXQ^^^ 

qSn+1 ~ V) "'^^ ^i^x = 5n + l) ^ for each x > 0. 

on^ (log nj^ on^ (log 77, 

— In the subcritical speed area (c > cq) there exists a constant K = K{c, f) > 0, such that, 
as n ^ 00, 

(lSn+1 ~ —1-T and P{Zx = 6n + l) r = ^ — - for each x > 0. 

77,"+^ Ac — Ac n'^ 



2 



Furthermore, qsn+k = P{Zx = 5n + /c) = for all n £ Z and k G {2, . . . ,6}. 

Remark 1.3. The idea of using singularity analysis for the study of Zx comes from Robin 
Pemantle's (unfinished) manuscript [24] about branching random walks with Bernoulli repro- 
duction. 

Remark 1.4. Since the coefficients of the power series i?[s'^] are real and non-negative, Pring- 
sheim's theorem (see e.g. [13], Theorem IV. 6, p. 240) entails that the assumption in Theorem 
1.2 is verified if and only if f{s) is analytic at 1. 

Remark 1.5. Let /? > and o" > 0. We consider a more general branching Brownian motion 
with branching rate given by /3 and the drift and variance of the Brownian motion given by 
c and (T^, respectively. Call this process the (/3, c, cr)-BBM (the reproduction is still governed 
by the law of L, which is fixed). In this terminology, the process described at the beginning 
of this section is the (1, c, 1)-BBM. The (/3, c, o-)-BBM can be obtained from {l,c/{ay/^), 1)- 
BBM by rescaling time by a factor /? and space by a factor cj/y^. Therefore, if we add 
an absorbing barrier at the point x > 0, the (/3, c, (t)-BBM gets extinct a.s. if and only if 

c > Co = a\/2f3m. Moreover, if we denote by zj^f'^''^^ the number of particles absorbed at x, 
we obtain that 

(Z(^'^''^)),>0 and {Z^^^^l'^^'''^)x>o are equal in law. 

In particular, if we denote the infinitesimal transition rates of {Zx^'^''^^)x>o by q^'^''^\ for 
n £ No\{l}, then we have 

One therefore easily checks that the statements of Theorems 1.1 and 1.2 are still valid for 
arbitrary [3 > and a > 0, provided that one replaces the constants cq, Xc, Xc, K by cq/ct^, 

Xc/a^ Va2, ^Kic/{a^),f), respectively. 

Remark 1.6. After submission of this article, Yang and Ren published an article [25] which 
permits to weaken the hypothesis in Theorem 1.1: It is enough to assume that -E'[L(log L)"^] < 
oo. In our proof, one needs to replace the reference [20] by [25] and use Theorem B of [7] 
instead of our Lemma 4.1, in order to obtain (4.6). 

The content of the paper is organised as follows: In Section 2 we derive some preliminary 
results by probabilistic means. In Section 3, we recall a known relation between Zx and the 
so-called Fisher-Kolmogorov-Petrovskii-Piskounov (FKPP) equation. Section 4 is devoted 
to the proof of Theorem 1.1, which draws on a Tauberian theorem and known asymptotics 
of travelling wave solutions to the FKPP equation. In Section 5 we review results about 
complex differential equations, singularity analysis of generating functions and continuous- 
time Galton-Watson processes. Those are needed for the proof of Theorem 1.2, which is done 
in Section 6. 

2 First results by probabilistic methods 

The goal of this section is to prove 
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Proposition 2.1. Assume c> cq and E[L'^] < oo. There exists a constant C = C{x,c,L) > 
0, such that 

C 

P{Zx > n) > — for large n. 

This result is needed to assure that the constant K in Theorem 1.2 is non-zero. It is 
independent from Sections 3 and 4 and in particular from Theorem 1.1. Its proof is entirely 
probabilistic and follows closely [2]. 

2.1 Notation and preliminary remarks 

Our notation borrows from [20]. An individual is an element in the space of Ulam-Harris 
labels 

U= J N", 

neNo 

which is endowed with the ordering relations ^ and -< defined by 

u <v <^=^ 3w £ U : V = uw and u ~< v u ^ v and u ^ v. 

The space of Galton-Watson trees is the space of subsets i C U, such that G t, -y G t if 
V -< u and u £ t and for every u there is a number G Nq, such that for all j G N, uj G t if 
and only if j < L^. Thus, is the number of children of the individual u. 

Branching Brownian motion is defined on the filtered probability space {T , J- , (J-t) , P) ■ 
Here, T is the space of Galton-Watson trees with each individual u G t having a mark 
{Cu,Xu) G M+ X i:'(M+,MU{A}), where A is a cemetery symbol and i:)(M+, MU {A}) denotes 
the Skorokhod space of cadlag functions from M"*" to MU {A}. Here, denotes the life length 
and Xu{t) the position of u at time t, or of its ancestor that was alive at time t. More precisely, 
for V £ t, let dy = J2w~<v denote the time of death and by = dy — the time of birth of v. 
Then Xu{t) = A for t > dy and if v ^ u is such that t G [by,dy), then Xu{t) = Xy{t). 
The sigma-field J-f contains all the information up to time t, and T = a ^Ut>o -^i) • 
Let y, c G M and L be some random variable taking values in No\{l}. P = PV''^'^ is the 
unique probability measure, such that, starting with a single individual at the point 

— each individual moves according to a Brownian motion with drift c until an independent 
time C,u following an exponential distribution with parameter 1. 

— At the time Cm the individual dies and leaves L„ offspring at the position where it has 
died, with being an independent copy of L. 

— Each child of u repeats this process, all independently of one another and of the past of 
the process. 

Note that often c and L are regarded as fixed and y as variable. In this case, the notation P^ 
is used. In the same way, expectation with respect to P is denoted by E or . 

A common technique in branching processes since [21] is to enhance the space T by 
selecting an infinite genealogical line of descent from the ancestor 0, called the spine. More 
precisely, if T G T and t its underlying Galton-Watson tree, then ^ = (^Oj Cij ^2, • • •) S is 
a spine oiT \i = and for every n G Nq, Cn.+i is a child of ^„ in t. This gives the space 

f = {(T, G r X C/^o : ^ is a spine of T] 

of marked trees with spine and the sigma-fields J- and Tt- Note that if (T, ^) G T, then T is 
necessarily infinite. 
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Assume from now on that m = E[L] — 1 G (0, oo). Let Nt be the set of individuals ahve 
at time t. Note that every /"^-measurable function / : T — )• M admits a representation 

where is an J-j-measurable function for every u £ U. We can therefore define a measure P 
on {f,T, {Ft)) by 

/ fdP = e~"'' [ Y.UT)P{dT). (2.1) 

It is known [20] that this definition is sound and that P is actually a probability measure 
with the following properties: 

— Under P, the individuals on the spine move according to Brownian motion with drift c 
and die at an accelerated rate m + 1, independent of the motion. 

— When an individual on the spine dies, it leaves a random number of offspring at the 
point where it has died, this number following the size-biased distribution of L. In other 
words, let L be a random variable with E[f{L)] = E[f{L)L/{m + 1)] for every positive 
measurable function /. Then the number of offspring is an independent copy of L. 

— Amongst those offspring, the next individual on the spine is chosen uniformly. This 
individual repeats the behaviour of its parent. 

— The other offspring initiate branching Brownian motions according to the law P. 
Seen as an equation rather than a definition, (2.1) also goes by the name of "many-to-one 
lemma" . 



2.2 Branching Brownian motion with two barriers 

We recall the notation P^ from the previous subsection for the law of branching Brownian 
motion started at y G M and the expectation with respect to P^ . Recall the definition of 
P and define P^ and E^ analogously. 

Let a, 6 G M such that y £ {a,b). Let r = Ta^b be the (random) set of those individuals 
whose paths enter (— oo, a] U [b, oo) and all of whose ancestors' paths have stayed inside (a, b). 
For u £ T we denote by r(u) the first exit time from (a, b) by u's path, i.e. 

r(u) = inf{i > : X„(t) ^ (a, 6)} = mm{t > : X„(t) G {a, 6}}, 

and set t{u) = oo for u ^ t. The random set r is an (optional) stopping line in the sense of 
[10]. 

For u £ T, define X„(t) = Xu{t{u)). Denote by Za^b the number of individuals leaving 
the interval (a, b) at the point a, i.e. 

Za,b = Y lXu{r)=a- 



Lemma 2.2. Assume |c| > cq and define p = \J (? — Cq . Then 

smh((o — a)p) 
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//, furthermore, V = E[L{L — 1)] < oo, then 
2ye'={«-s/) 



psinh ((5 — a)p) 



ry 

sinh((6 - / e''^""''^ sinh^((6 - r)p) smh((r - a)p) dr 

J a 



+ smh((y-a)p) / e''^''"''hinh\{b - r)p) dr] + Ey[Za,b]- 



Proof. On the space T of marked trees with spine, define the random variable / by / = i if 
£,i £ T and I = oo otherwise. For an event A and a random variable Y write E\Y, A] instead 
ofE[YlA]. Then 

Ey[Za,t] = Ey[Y,^x^ir)=a] = [e™'^(«^\ / < oo,X^^{t) = a] 

u£t 

by the many-to-one lemma extended to optional stopping lines (see [6], Lemma 14.1 for a 
discrete version). But since the spine follows Brownian motion with drift c, we have / < oo, 
P-a.s. and the above quantity is therefore equal to 

where W^''^ is the law of standard Brownian motion with drift c started at y, {Bt)t>o the 
canonical process and T = Ta^b the first exit time from (a, 6) of Bt. By Girsanov's theorem, 
and recalling that m = this is equal to 

where = W^'^ . Evaluating this expression ([8], p. 212, Formula 1.3.0.5) gives the first 
equality. 

For u £ U, let 0„ be the operator that maps a tree in T to its sub-tree rooted in u. 
Denote further by C„ the set of tt's children, i.e. C„ = {uk : 1 < k < L^}. Then note that for 
each u £ T we have 

Za,b = 1 + ^a,b ° &w, 

hence 



Ey[zl,] = Ey[Y^ix^f^,^=aZa,i 



EyiZaA+P" 



(2.2) 



Define the a- algebras 



n = gva{c^;v -<Ci), 
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such that Q contains the information about the path of the spine up to the individual that 
quits (a, b) first, H adds to G the information about the fission times on the spine and I adds 
to Ti the information about the individuals of the spine and the number of their children. 
Now, conditioning on I and using the strong branching property, the second term in the last 
line of (2.2) is equal to 



py 



(recall that is the time of death of v). Conditioning on % and noting the fact that L„ 
follows the size-biased law of L for an individual v on the spine, yields 



py 



V 



Finally, since under P the fission times on the spine form a Poisson process of intensity m + 1, 
conditioning on G and applying Girsanov's theorem yields 



^c{BT~y)~\p^T 



r VE^'^lZa^b] dt,BT 
Jo 



dr. 



where is the local time of (Bt) at the time T and the point r. The last expression can be 
evaluated explicitly ([8], p. 215, Formula 1.3.3.8) and gives the desired equality. □ 

Corollary 2.3. Under the assumptions of Lemma 2.2, for each b > there are positive 
constants cj^^\ such that as a ^ —oo, 

a) i?0[Z,,,]~C^')e(^+^)«, 

b) ifo Co, E^[Zlf;\ ~ cl^^e^^+P^" and 

c) ifc<-co, e2(-+'')^ 

The following result is well known and is only included for completeness. We emphasize 
that the only moment assumption here is m = E[L] — 1 G (0,oo). Recall that Z^ denotes the 
number of particles absorbed at x of a BBM started at the origin. For |c| > cq, define Ac to 
be the smaller root of — 2c + Cq , thus Ac = c — y^c^ ~ Cg . 



XcX 



Lemma 2.4. Let x > 0. 

— If |c| > Co, then E[Zx] = e 

— // |c| < Co, then E[Zx] = +oo. 

Proof. We proceed similarly to the first part of Lemma 2.2. Define the (optional) stopping 
line r of the individuals whose paths enter [x, oo) and all of whose ancestors' paths have 
stayed inside (—00,2;). Define / as in the proof of Lemma 2.2. By the stopping line version 
of the many-to-one lemma we have 

E[Zx]=E[Y^ l] = E[e'"^«^),/<oo]. 
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By Girsanov's theorem, this equals 



Ty[e=^-^('='-'=o)^-,T, < oo], 

where W is the law of standard Brownian motion started at and is the first hitting time 
of X. The result now follows from [8], p. 198, Formula 1.2.0.1. □ 

2.3 Proof of Proposition 2.1 

By hypothesis, c > cq, E[L'^] < oo and the BBM starts at the origin. Let x > and let t = Tx 
be the stopping line of those individuals hitting the point x for the first time. Then Zx = \tx\- 
Let a < and n G N. By the strong branching property, 

P\Zx >n)> P\Zx > n I Za,x > l)P°(^a,x- > 1) > P'^iZx > n)P\Za,x > 1). 

If denotes the law of branching Brownian motion started at the point with drift — c, 
then 

P^{Zx >n) = Pl{Za-x >n)> Pl{Za-x,i > n). 

In order to bound this quantity, we choose a = a„ in such a way that 7i = ^E^[Za„-x,i]- By 
Corollary 2.3 a), c) (applied with drift — c) and the Paley-Zygmund inequality, there is then 
a constant Ci > 0, such that 

P-{Zar,-x,i > n) > j^^f^^^ > Ci for large n. 

Furthermore, by Corollary 2.3 a) (applied with drift — c), we have 

^Cf ^e-^=(""-^) ~ n, as n ^ oo, 

and therefore an = — (1/Ac) log n + 0(1). Again by the Paley-Zygmund inequality and Corol- 
lary 2.3 a), b) (applied with drift c), there exists C2 > 0, such that for large n, 

P (Z.„, > 1) - P (Z.„, > 0) > > — ^e > 

This proves the proposition with C = C1C2. 

3 The FKPP equation 

As was already observed by Neveu [23], the translational invariance of Brownian motion 
and the strong branching property immediately imply that Z = {Zx)x>o is a homogeneous 
continuous-time Galton- Watson process (for an overview to these processes, see [4], Chap- 
ter III or [14], Chapter V). There is therefore an infinitesimal generating function 

a{s) = ay^PnS^ - s\ , Q>0, pi=0, (3.1) 

\n=0 / 
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associated to it. It is a strictly convex function on [0,1], with a(0) > and a(l) < 0. Its 
probabilistic interpretation is 

a = lim \P{Zr^ / 1) and p„ = lim P{Z^ = n\Z^ / 1), 

x^O x^O 

hence g„ = apn for n G No\{l}. Note that with no further conditions on c and L, the sum 
X]n>o^'" need not necessarily be 1, i.e. the rate apoo, where poo = 1 — Z^,i>oPni with which 
the process jumps to +oo, may be positive. 

We further define Fx{s) = Els^"^], which is linked to a{s) by Kolmogorov's forward and 
backward equations ([4], p. 106 or [14], p. 102): 

d d 

—Fx(s) = a(s) — Fx(s) (forward equation) (3.2) 

ox OS 

d 

—Fx(s) = a[Fx(s)] (backward equation) (3.3) 

ox 

The forward equation implies that if a(l) = and (j){x) = E[Zx] = -^Fx{l — ), then (j)'{x) = 
a'{l)(j){x), whence E[Zx] = e'^'(^)^. On the other hand, if a(l) < 0, then the process jumps to 
oo with positive rate, hence E[Zx] = oo for all x > 0. 

The next lemma is an extension of a result which is stated, but not proven, in [23], 
Equation (1.1). According to Neveu, it is due to A. Joffe. To the knowledge of the author, 
no proof of this result exists in the current literature, which is why we prove it here. 

Lemma 3.1. Let (lt)j>o be a homogeneous Galton-Watson process started at 1, which may 
explode and may jump to +oo with positive rate. Let u{s) be its infinitesimal generating 
function and Ft{s) = ii^[s^*]. Let q be the smallest zero of u{s) in [0, 1]. 

1. If q < 1, then there exists t_ G M U {— oo} and a strictly decreasing smooth function 
ifj- : (t_,+oo) — {q,l) with limt^t^ fp-{t) = 1 and limt^oo i^-{t) = q, such that on 
{q, 1) we have u = ip'_ o ipZ^ , Ft{s) = ip-{'ipZ^{s) + t). 

2. If q > 0, then there exists t^ G ML) { — oo} and a strictly increasing smooth function 
V'+ : (t+,+oo) — (0, g) with limt-^t+ ip+{t) = and limt^oo i^+it) = Q, such that on 
(0, q) we have u = o i})'^^ , Ft{s) = '(/'+('0+"'"(s) + t). 

The functions and ipj^ are unique up to translation. 
Moreover, the following statements are equivalent: 

— For all t > 0, Yf < oo a.s. 

— q = 1 or = — oo. 

Proof. We first note that u{s) > on {0,q) and u{s) < on ((7,1), since u{s) is strictly 
convex, ii(0) > and n(l) < 0. Since Fq{s) = s, Kolmogorov's forward equation (3.2) implies 
that Ft{s) is strictly increasing in t for s G (0, (?) and strictly decreasing in i for s G ((?, 1)- 
The backward equation (3.3) implies that Ft{s) converges to (? as t — )■ oo for every s G [0, 1). 
Repeated application of (3.3) yields that Ft{s) is a smooth function of t for every s G [0, 1]. 

Now assume that q < 1. For n G N set s„ = 1 — 2~"(1 — q), such that q < si < 1, 
Sn < Sn+i and s„ — )• 1 as ?i — )• oo. Set ti = and define t„ recursively by 

tn+i = tn — i', where t' > is such that Fti{sn+i) = Sn- 

Then {tn) ngN is a decreasing sequence and thus has a limit t_ G M U {— oo}. We now define 
for t G (t_, +oo), 

ij^{t) = Ft^t,Xsn), ift>in. 
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The function ip- is well defined, sincG for every ?i G M ajiid t ^tni 

by the branching property. The same argument shows us that if s G {q, 1), s„ > s and t' > 
such that Fti{sn) = s, then Ft{s) = Ft^t'{sn) = i^-it + t' + tn) for all t > 0. In particular, 
ip-{t' + tn) = s, hence Ft{s) = i{j^{'ipZ^{s) + 1). The backward equation (3.3) now gives 

u{s) = ^^Ft{s)^^^^ = ^P'_{^PZ\s)). 

The second part concerning tp^ is proven completely analogously. Uniqueness up to transla- 
tion of tjj- and ■0+ is obvious from the requirement ip{tp~^{s) + 1) = Ft{s), where tp is either 
tp- or ip^. 

For the last statement, note that P(Yt < oo) = 1 for all t > if and only if Ft{l—) = 1 
for all t > 0. But this is the case exactly if g = 1 or t_ = — oo. □ 

The following proposition shows that the functions ip- and ip+ corresponding to {Zx)x>o 
are so-called travelling wave solutions of a reaction-diffusion equation called the Fisher- 
Kolmogorov-Petrovskii-Piskounov (FKPP) equation. This should not be regarded as a new 
result, since Neveu ([23], Proposition 3) proved it already for the case c > cq and L = 2 a.s. 
(dyadic branching). However, his proof relied on a path decomposition result for Brownian 
motion, whereas we show that it follows from simple renewal argument valid for branching 
diffusions in general. 

Recall that /(s) = -E[,s^] denotes the generating function of L. Let q' be the unique fixed 
point of / in [0, 1) (which exists, since /'(I) = m + 1 > 1), and let q be the smallest zero of 
a{s) in [0,1]. 

Proposition 3.2. Assume c G R. The functions ■0- and 0+ from Lemma 3.1 corresponding 
to {Zx)x>o o,re solutions to the following differential equation on (t_,+oo) and (t_|_,+oo), 
respectively. 

h^" -ci;' = i;-fo^. (3.4) 

Moreover, we have the following three cases: 

1. If c> Co, then q = q' , t_ = -oo, a(l) = 0, a'(l) = Ac, E[Zx] = e'^'^^ for all x > 0. 

2. If \c\ < Co, then q = q' , t„ G M, a(l) < 0, a'(l) = 2c, P{Zx = oo) > for all x > 0. 

3. Ifc< -Co, then q = 1, a(l) = 0, a'(l) = Ac, E[Zx] = e^'^^ for all x > 0. 

Proof. Let s G (0,1) and define the function ips{x) = Fx{s) = £'[s^^] for x > 0. By sym- 
metry, Zx has the same law as the number of individuals absorbed at the origin in a 
branching Brownian motion started at x and with drift — c. By a standard renewal argument 
(Lemma A.l), the function ipg is therefore a solution of (3.4) on (0, oo) with ■0s(O+) = s. This 
proves the first statement, in view of the representation of Fx in terms of ^— and ^p^ given 
by Lemma 3.1. 

Let s G (0, 1)\{(?} and let ip{s) = tp-is) if s > q and tpis) = tp+{s) otherwise. By (3.4), 

, ip"oip-^{s) %p o^-^{s) - f oip oip-''{s) ,^s-f{s) 
« Os) = — — =2c + 2 — — -— = 2c + 2 -— , 

ip'oip ^[s) ip'oip ^[s) a{sj 
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whence, by convexity, 

a'{s)a{s) = 2ca{s) + 2{s- f{s)), sG[0, 1]. (3.5) 

Assume |c| > cq. By Lemma 2.4, E[Zx] = e'^''^, hence o(l) = and a'(l) = Ac, in 
particular, a'(l) > for c > cq and a'(l) < for c < —cq. By convexity, g < 1 for c > cq and 
q = 1 for c < — Cq. The last statement of Lemma 3.1 now implies that t_ = — oo if c > cq. 

Now assume |c| < cq. By Lemma 2.4, E[Zx] = +oo for all a; > 0, hence either a(l) < 
or a(l) = and a'(l) = +oo, in particular, g < 1 by convexity. However, if a(l) = 0, then 
by (3.5), a'(l) = 2c — 2m/a'(l), whence the second case cannot occur. Thus, a(l) < and 
a'(l) = 2c by (3.5). 

It remains to show that q = q' if q < 1. Assume q ^ q' . Then a{q') 7^ by the (strict) 
convexity of a and a'{q') = 2c by (3.5). In particular, a'{q') > a'{l), which is a contradiction 
to a being strictly convex. □ 



4 Proof of Theorem 1.1 

We start with the following Abelian-type lemma: 

Lemma 4.1. Let X be a random variable concentrated on No and let (/^(s) = i?[s'''"] be its 
generating function. Assume that E[X(log~^ -^V] < oo for some 7 > 0. Then, as s — )■ 0, 

<^'(l)_^'(l_5) = 0((logl)-7) and ^'{l)s + ^{l-s)-l = Ois{log^r^). 

Proof. Let sq > be such that the function s 1— )• s(logi)''' is increasing on [0, sq]. Let 
s £ (0, So)- Then, with pk = P{X = k), 

00 

((^'(1) - ^'(1 - s))(iog i)^ = kPkil - (1 - s)'-^)(iog ir. 

k=l 

If fe > s-\ then (1 - (1 - s)'=-i)(logi)^ < (logA;)^. If [sq^] < A: < s'^, then s(logi)T < 
lilogk)'^ and thus (1 - (1 - sf~^){log^)'^ < A:s(logi)T < (logA;)'^. Hence, 

00 00 

A:pfc(l-(l-s)'"^)(logi)^< Y PkH^ogky <E[X{log+Xn 

k=\s-^] k=\s-^] 

Furthermore, we have for s G (0, 1), 

J2 kpkii - (1 - s)'-')(iog '-r < (s(iog i)^) Y ^v. < c, 

fc=l fc=l 
for some C > 0. Collecting these results, we have, for every s G (0, 1), 

(<^'(1) - ip'il - s))(log i)^ < C + i?[X(log+ Xy] < oo, 

by hypothesis. This yields the first equality. Setting g{s) = ^'{l)s + ip{l — s) — 1, we note 
that g{0) = and 

5'(s) = ^'(l)-V^'(l-s) = 0((logl)"^), 
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by the first equality. Since (log ^ ) ^ is slowly varying, 

gis)= r<7'(r)dr = 0(s(logl)-^) 



by standard theorems on the integration of slowly varying functions (see e.g. [11], Section 
VIII.9, Theorem 1). □ 

Proof of Theorem 1.1. We have c = cq by hypothesis. Let ip- be the travelling wave from 
Proposition 3.2, which is defined on M, since t_ = — oo. Let (/>(x) = 1 — ^/;_(— x), such that 
(^(— oo) = 1 — g, <^(+oo) = and 

y\x) + cocp'ix) = /(I - 0(x)) - (1 - cPix)), (4.1) 

by (3.4). Furthermore, a(l — s) = <p' {(f)^^ [s]) and Fx{l — s) = 1 — 0((/>~^(s) — x). 

Under the hypothesis £'[L(log L)^+"] < oo, it is known [20] that there exists K e (0, oo), 
such that (j){x) ~ Kxe~'^°^ as x ^ oo. Since a(l) = and a'(l) = cq by Proposition 3.2, this 
entails that (f>'{x) = a(l — (j){x)) ~ —coKxe~'^°^ , as x — )■ oo. 

Set = (j)' and (p2 = 4>. By (4.1), 

A = f'i'"i^)\ = ^-2co0'(x) + 2[/(l - m) - (1 - cPix))] 

dx\<f2{x)J U'{x)J V <t>'ix) 

Setting g{.s) = cgs + 2[/(l - s) - (1 - s)] = 2[f (l)s + /(I - s) - 1], this gives 



-giHx)) 



dx \(f2{x) J \(f2{x) J \ J' \ 1 

The Jordan decomposition of M is given by 

J-^-MA.{7 -co)' ^-("" "l")- 
Setting (^;)=^(|)'-<^g«twithe=(|): 

e'(x) = je(x) + 

which, in integrated form, becomes 

«.) =e-«0) + e-£e-.^ d,. (4.4) 

Note that 

By the definition of ^ and the above asymptotic of (p, we have g{(j){x)) = 0{e~'^°^ /x^'^'^), 
as X — )• 00, by Lemma 4.1 and the hypothesis on L. Equations (4.4) and (4.5) now imply that 

6(x) ~ e-'=o^(6(0) + r e-'yg{<P{y))dy), 
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and 

(6 + ^2){x) ~ xe-=°^ (6(0) + e=«^'5(</'(y)) dy) , 
and since = + 6 and 6.2 = (p' + CQfj), this gives 

{(p' + co^){x) (t){x)/x Ke-''°''. (4.6) 
With this information, one can now show by elementary calculus (see Section A. 2), that 



- ~ a^^d (4.7) 



_co 
s(log 
coxe""^"^ 



i^'(l - s) ~ ^^^f^, as . ^ 0. (4.8) 

S(logi)^ 

By standard Tauberian theorems ([11], Section XIII.5, Theorem 5), (4.7) implies that 

n 

U{n) = y^ k^qk ^ cq- -k, as n 00. 

^ (logn)^ 

By integration by parts, this entails that 

^^'^ Jn-^ f^(da;) x2(logx)2'^^ re(logn)2j " 

But the last integral is equivalent to l/(n(logn)^) ([11], Section VIII. 9, Theorem 1), which 
proves the first part of the theorem. The second part is proven analogously, using (4.8) 
instead. □ 

5 Preliminaries for the proof of Theorem 1.2 

In light of Proposition 2.1, one may suggest that under suitable conditions on L one may 
extend the proof of Theorem 1.1 to the subcritical case c > cq and prove that as n — )• 00, 
P{Zx > n) ~ C'n~'^ for some constant C". In order to apply Tauberian theorems, one would 
then have to establish asymptotics for the {[d\ + l)-th derivatives of a(s) and Fx{s) as s — )• 1. 
In trying to do this, one quickly sees that the known asymptotics for the travelling wave 
(1 — Tp{x) ~ const X e""^"^ as x — )• —00, see [20]) are not precise enough for this method to 
work. However, instead of relying on Tauberian theorems, one can analyse the behaviour 
of the holomorphic function a{s) near its singular point 1. This method is widely used in 
combinatorics at least since the seminal paper by Flajolet and Odlyzko [12] and is the basis 
for our proof of Theorem 1.2. Not only does it work in both the critical and subcritical cases, 
it even yields asymptotics for the density instead of the tail only. 

In the rest of this section, we will define our notation for the complex analytic part of the 
proof and review some necessary general complex analytic results. 

5.1 Notation 

In the course of the paper, we will work in the spaces C and C^, endowed with the Euclidean 
topology. An open connected set is called a region, a simply connected region containing a 
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point zq is also called a neighbourhood of zq. The closure of a set D is denoted by D, its 
border by dD. The disk of radius r around zq is denoted by D(zo, r) = {z £ C : \z — zo\ < r}, 
its closure and border by ^{zq, r) and dD{zo,r), respectively. We further use the abbreviation 
D = D(0, 1) for the unit disk. For 0<93<7r, r>0 and x G M, we define 

G{ip, r) = {z G D(l, r)\{l} : | arg(l - 2;)| < vr - c^}, S+{ip, x) = [x, oo) x {-cp, ip), 

A(v9,r) = {ze D(0,1 + r)\{l} : | arg(l - z)\ < vr - (^}, S^{ip,x) = (-00, x] x {-if, if), 
H{ip, r) = {z G D(0, r)\{0} : | arg z\ < ip}. 

Note that H{{p, r) = 1 — G{tt — 99, r). Here and during the rest of the paper, arg(z) and log(z) 
are the principal values of argument and logarithm, respectively. 

Let G be a region in C, zq £ G and / and g analytic functions in G with g{z) 7^ for all 
z £ G. We write 

f{z) = o{g{z)) ^ ye> 36 >OyzeGnO{zo,6):\f{z)\<£\g{z)\, 
f{z) = 0{g{z)) ^ 3C >0 35 > Vz G GnD(zo,<^) : 1/(^)1 <C|5(2)|, 
f{z) = d{g{z)) ^ 3K€C: f{z) = Kg{z) + o{g{z)), 
f{z)^g{z) ^ f{z)=g{z) + o{g{z)), 

specifying that the relations hold as 2; — >■ zq. 



5.2 Complex differential equations 

In this section, we review some basics about complex differential equations. We start with 
the fundamental existence and uniqueness theorem ([5], p. 1, [15], Theorem 2.2.1, p. 45 or 
[18], Section 12.1, p. 281). 

Fact 5.1. Let G he a region in and (wq, zq) a point in G. Let f : G ^ C be analytic in G, 
i. e. f is continuous and both partial derivatives exist and are continuous. Then there exists a 
neighbourhood U of zq and a unique analytic function w : U ^ C, such that 

1. w{zq) = Wo, 

2. {w{z),z) G G for all z £U and 

3. w'{z) = f{w{z),z) for all z £U. 

In other words, the differential equation w' = f{w,z) with initial condition w{zo) = wq has 
exactly one solution w{z) which is analytic at zq. 

The following standard result is a special case of a theorem by Painleve ([5], p. 11, [15], 
Theorem 3.2.1, p. 82 or [18], Section 12.3, p. 286f). 

Fact 5.2. Let H be a region in C and w{z) analytic in H. Let G be a region in C^, such that 
{w{z),z) G G for each z £ H and suppose that there exists an analytic function / : G — >■ C, 
such that w'{z) = f{w{z),z) for each z £ H. Let zq G dH . Suppose that w{z) is continuous 
at Zq and that {w{zo), zq) G G. Then zq is a regular point ofw{z), i.e. w{z) admits an analytic 
extension at zq. 

Let [ 

^1) • • • ) ^k]n denote a power series of the variables zi, . . . , z^, converging in a neigh- 
bourhood of (0, . . . , 0) and which contains only terms of order n or higher. The complex 
differential equation 

zw' = Xw + pz + [w, z]2, A,p G C, (5-1) 
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was introduced in 1856 by Briot and Bouquet [9] as an example of a complex differential 
equation admitting analytic solutions at a singular point of the equation. More precisely, 
they obtained ([15], Theorem 11.1.1, p. 402): 

Fact 5.3. // A is not a positive integer, then there exists a unique function w{z) which is 
analytic in a neighbourhood of z = and which satisfies (5.1). Furthermore, w{0) = 0. 

The singular solutions to this equation were later investigated by Poincare, Picard and 
others (for a full bibliography, see [17]). We are going to need the following result (see [17], 
Paragraph 111.9.2° or [15], Theorem 11.1.3, p. 405, but note that the latter reference is without 
proof and the statement is slightly incomplete). 

Fact 5.4. Assume A > 0. There exists a function ijj{z,u) = ^jk>oPjkZ-'u^ , converging in a 
neighbourhood of (0,0) and such that poo = and poi = 1? such that the general solution of 
(5.1) which vanishes at the origin is w = 'ip{z,u), with 

- u = Cz^, ifXin, 

- u = z^{C + K\ogz), ifXe N. 

Here, C ^ C is an arbitrary constant and K & C is a fixed constant depending only on the 
right-hand side of (5.1). 

Remark 5.5. The above statement is slightly imprecise, in that the term solution is not 
defined, i.e. what a priori knowledge of w{z) (regarding its domain of analyticity, smoothness, 
behaviour at z = 0, . . . ) is required in order to guarantee that it admits the representation 
stated in Fact 5.4? Inspecting the proof (as in [17], for example) shows that it is actually 
enough to know that w{z) satisfies (5.1) on an interval (0, e) of the real line and that w{0+) = 
0. We briefly explain why: 

In order to prove Fact 5.4, one shows that there exists a function of the form stated 
above, such that when changing variables hy w = ip{z, u), the function u{z) formally satisfies 
one of the equations 

zu' = An or zu' = Xu + Kz^, 

according to whether A ^ N or A G N. 

Now suppose that w{z) satisfies the above conditions. By the implicit function theorem 
([16], Theorem 2.1.2), we can invert V to obtain a function Lp{w,z) = w + qz + [w,z\2-, 
q £ C, such that ip{z,ip{w, z)) = in a neighbourhood of (0,0). We may thus define 
u{z) = ip{w{z),z) for all z € (0,ei) for some ei > 0. Moreover, u{z) now truly satisfies the 
above equations on (0,ei) and n(0+) = 0. Standard theory of ordinary differential equations 
on the real line now yields that u is necessarily of the form stated in Fact 5.4. 

We further remark that since u{z) is analytic in the slit plane C\(— oo, 0] and goes to as 
z — )• in C\(— oo, 0], there exists an r > 0, such that (z, u{z)) is in the domain of convergence 
of ip{z,u) for every z £ H{Tr,r). Hence, every solution w{z) can be analytically extended to 
H{7T,r). 

5.3 Singularity analysis 

We now summarise results about the singularity analysis of generating functions. The basic 
references are [12] and [13], Chapter VI. The results are of two types: those that establish an 
asymptotic for the coefficients of functions that are explicitly known, and those that estimate 
the coefficients of functions which are dominated by another function. We start with the 
results of the first type: 
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Fact 5.6. Let d G (l,oo)\N, /c G N, 7 G Z\{0}, 5 G Z and the functions fi, f2 defined by 



h{z) = il-zy and /2(z) = (l-2)'=(^log^y (^loglog^) , 

for z G C\[l,+oo). Let {pn^) be the coefficients of the Taylor expansion of fi around the 
origin, i = 1,2. Then {p^n) satisfy the following asymptotics as n ^ 00; 



(1) K, i^2(logn)^-Hloglogn)'^ 



for some non-zero constants Ki = Ki{d),K2 = K2{k,j,5). We have K2{1, —1,0) = 1. 

Proof. For /i, this is Proposition 1 from [12]. For /2 this is Remark 3 at the end of Chapter 
3 in the same paper. Note that the additional factors ^ do not change the nature of the 
singularities, since ^ is analytic at 1 (see the footnote on p. 385 in [13]). The last statement 
follows from Remark 3 as well. □ 

The results of the second type are contained in the next theorem. It is identical to 
Corollary 4 in [12]. Note that a potential difficulty here is that it requires analytical extension 
outside the unit disk. 

Fact 5.7. Let < ip < tt/2, r > and f{z) be analytic in A{Lp,r). Assume that as z — )• 1 in 
A{ip,r), 

/(z) = 0(^(1 -z)"L^^-^^^ , w/iere L('u) = (logn)^(loglogn)'', a,j,S eR. 
Then the coefficients (p„,) of the Taylor expansion of f around satisfy 

Pn = o\ — — r , as CO. 



5.4 An equation for continuous-time Galton— Watson processes 

In this section, let (lt)f>o be a homogeneous continuous-time Galton- Watson process starting 
at 1. Let a{s) be its infinitesimal generating function and Ft{s) = E[s^^]. Assume a(l) = 
and a'(l) = A G (0, 00), such that a(s) = has a unique root q in [0, 1). 

The following proposition establishes a relation between the infinitesimal generating func- 
tion of a Galton-Watson process and its generating function at time t. For real s, the formulae 
stated in the proposition are well known, but we will need to use them for complex s, which 
is why we have to include some (complicated) hypotheses to be sure that the functions and 
integrals appearing in the formulae are well defined. 

Proposition 5.8. Suppose that a and Ft have analytic extensions to some regions Da and 
Dp. Let Za = {s £ Da : a{s) = 0}. Let there be simply connected regions G C Da\Za and 
D C G D Dp with Ft{D) C G and D H (0, 1) 7^ 0. Then the following equations hold for all 
s G D: 

'■Ft{s) I 

^dr = t, (5.2) 



s 
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and 



l-Ft{s)=e^\l-s)expl- J r(r)drj 



(5.3) 



where f*{s) is defined for all s G Da\Za as 



r(s) = ^ + :r^, (5.4) 
a[sj 1 — s 

and the integrals may be evaluated along any path from s to Ft{s) in G. 

Proof. For s G (0, l)\{g}, equation (5.2) follows readily from Kolmogorov's backward equation 
(3.3), when the integral is interpreted as the usual Riemann integral ([4], p. 106). Now note 
that by definition of G, both and /* are analytic in the simply connected region G and 
therefore possess antiderivatives g and h in G. Thus, the functions 



s H> 



/ dr = g{Ft{s)) - gis) and f*{r) dr = h{Ft{s)) - h{s) 



are analytic in D. By the analytic continuation principle, (5.2) then holds for every s € D, 
since -D n (0, 1) ^ by hypothesis. This proves the first equation. For the second equation, 
note that — log(l — s) is an antiderivative of in G, whence the right-hand side of (5.3) 
equals 



e^*(l - s) exp ^log(l - Ft{s)) - log(l " " ^ ^ 



Ft is) 1 \ 

dr = 1 - F(s) 



a{r) 

for all s £D,hj (5.2). This gives (5.3). □ 

Corollary 5.9. If 1 is a regular point of a{s), then it is a regular point for Ft{s) for every 
t > 0. 

Proof. Define G = {s € D : Re s > q}. Then G D Za = 0, since q is the only zero of a in D 
(every probability generating function g with ^'(l) > 1 has exactly one fixed point q in D; 
this can easily be seen by applying Schwarz's lemma to o g o t, where r is the Mobius 
transformation of the unit disk that maps to g). Let si € {q, 1) be such that Ft{s) G G for 
every s £ H = {s £D : Res > si}. We can then apply Proposition 5.8 to conclude that (5.3) 
holds for every s G H. 

Since a{s) is analytic in a neighbourhood U of 1 by hypothesis, it is easy to show that /* 
is analytic in U as well. Thus, /* has an antiderivative F* m H L) U. We define the function 
g{s) = (1 — s)exp(F*(s)) on HUU. Since ^'(1) = — exp(F*(l)) ^ 0, there exists an inverse 
g~^ of g' in a neighbourhood Ui of ^(1) = 0. Let U2 C U he a^ neighbourhood of 1, such that 
e'*'*g(s) G Ui for every s G U2- Define the analytic function Ft{s) = g~^{e'^^g{s)) for s G 172- 
Then by (5.3), we have Ft{s) = Ft{s) for every s £ H Ci U2, hence Ft is an analytic extension 
of Ft at 1. □ 

Corollary 5.10. Suppose that a{s) has an analytic extension to G{(pQ,rQ) for some < 
(fQ < TT and ro > 0. Suppose further that there exist c G M, 7 > 1, such that a(l — s) = 
—As + Acs/ logs + 0{s/\ log si''') as s — )■ 0. Then for every Lpo < < n there exists r > 0, 
such that Ft{s) can be analytically extended to G{ip,r), mapping G{ip,r) into G(y?o,ro). 
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Proof. Recall that A > 0. By hypothesis, we can then assume that a(s) 7^ in G{<po,rQ) by 
choosing rg small enough. Then A/a has an antiderivative A on G{ipo,ro). Define B(s) = 
A{1 — s) for s G H{7r — ipQ,ro), such that 

„ = ^ + -^ + 0^ ' 



S(l - C/ logs + 0(1 log S|-T)) S slogs V^l logs|~™i'^(T'2) 

We can therefore apply Lemma A. 7 to B and deduce that there exist 931 G {ipQ,ip) and 
ri,r £ (0,ro), such that A is injective on G{ipi,ri) and such that A{s) + At G A{G{ipi,ri)) 
for every s G G{ip,r). Hence, -Ff(s) = A~^{A{s) + At) is defined and analytic on G{ip,r). By 
(5.2), Ft{s) = Ft{s) on G{ip,r) nD, hence -Ft is an analytic extension of Ft, mapping G{ip,r) 
into G(99i,ri) C G((/90)^o) by definition. □ 



6 Proof of Theorem 1.2 

We turn back to branching Brownian motion and to our Galton-Watson process Z = {Zx)x>o 
of the number of individuals absorbed at the point x. Throughout this section, we place 
ourselves under the hypotheses of Theorem 1.2, i.e. we assume that c > cq = \/2m and that 
the radius of convergence of /(s) = E[s^] is greater than 1. The equation A^ — 2cA + Cq = 
then has the solutions Ac = c — ^Jc^ — Cq and Ac = c + \J (? — Cq, hence Ac = Ac = cq if c = cq 
and Ac < Co < Ac otherwise. The ratio d = Xc/^c is therefore greater than or equal to one, 
according to whether c > cq or c = cq, respectively. Recall further that 5 G N denotes the 
span of L — 1. 

Let a{s) = Oi{^j^yQPkS^ — s) be the infinitesimal generating function of Z and let Fx{s) = 
Els^'']. We recall the equation (3.5) from Section 3: For s G [0, 1], 

o'(s)a(s) = 2ca(s) + 2(s - /(s)). (6.1) 

By the analytic continuation principle, this equation is satisfied on the domain of analyticity 
of a{s), in particular, on D. 

We now give a quick overview of the proof. Starting point is the equation (6.1). We are 
going to see that this equation is closely related to the Briot-Bouquet equation (5.1) with 
X = d. The representation of the solution to this equation given by Fact 5.4 will therefore 
enable us to derive asymptotics for a{s) near its singular point s = 1 (Theorem 6.4). Via the 
results in Section 5.4, we will be able to transfer these to the functions Fx{s) (Corollary 6.6). 
Finally, the theorems of Flajolet and Odlyzko in Section 5.3 yield the asymptotics for g„ and 
P{Zx = n). 

More specifically, we will see that the main singular term in the expansion of a(l — s) 
or Fx{l — s) near s = is s*^, if d ^ N and s'^logs, if d G N. At first sight, this dichotomy 
might seem strange, but it becomes evident if one remembers that we expect the coefficients 
of Fx{s) (i.e. the probabilities P{Zx = n), assume S = 1) to behave like l/n'^'^^, if d > 1 (see 
Proposition 2.1). In light of Fact 5.6, a logarithmic factor must therefore appear if d is a 
natural number, otherwise Fx{s) would be analytic at 1, in which case its coefficients would 
decrease at least exponentially. 

We start by determining the singular points of a(s) and Fx{s) on the boundary of the unit 
disk, which is the content of the next three lemmas. 

Lemma 6.1. Let X be a random variable with law {pk)keNo o.'i^d let 2; > 0. Then the spans 
of X — 1 and of Zx — 1 are equal to 5. 
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Proof. This follows from the fact that the BBM starts with one individual and the number 
of individuals increases by Z — 1 when an individual gives birth to / children. □ 



Lemma 6.2. If 6 = 1, then a{s) and {Fx{s))x>o are analytic at every sq £ dD\{l}. If S > 2, 
then there exist a function h{s) and a family of functions (Zix(s))x>0; analytic on D, such 
that 

a{s) = sh{s^) and Fx{s) = shx{s^), 
for every s G D. Furthermore, h and {hx)x>o o'^e analytic at every sq £ 5D\{1}. 
Proof. Assume first that 5 >2. Define 

h{s) = aC^Pi+Sns''-l) and = = 1 + 

n n 

By Lemma 6.1, Pk+Sn = P{Zx = k + 5n) = for every k € {2, . . . , 5} and n G Z, whence 
a{s) = sh{s^) and Fx{s) = shx{s^) for every s G D. 

We now claim that a and Fx are analytic at every sq G dB> with Sq ^ 1. Note that if 
S >2, this implies that h and hx are analytic at every sq G 5D\{1}, since the function s ^ 
has an analytic inverse in a neighbourhood of any s ^ 0. 

First note that by [11], Lemma XV. 2. 3, p. 475, we have | J2nPn^o\ ^ ^ every sq G dO, 
such that Sq 7^ 1, whence a{so) ^ 0. Now write the differential equation (6.1) in the form 

/ 2ca + 2{s-f{s)) 

a = =: g{a, s). 

a 

Since the radius of convergence of / is greater than 1 by hypothesis, g is analytic at (a(so), so)- 
Furthermore, a is continuous at sq, since "^nPnS^ converges absolutely for every s G O. 
Fact 5.2 now shows that a is analytic at sq. 

It remains to show that Fx is analytic at sq. Kolmogorov's forward and backward equations 
(3.2) and (3.3) imply that a(s)F^(s) = a{Fx{s)) on [0,1], and the analytic continuation 
principle implies that this holds on D. Now, let sq G (9D, such that Sq 1. Then we have 
just shown that a is analytic and non-zero at sq. Furthermore, |-Fa;('So)| < 1) by the above 
stated lemma in [11] and Lemma 6.1. Thus, the function f{w,s) = a{w)/a{s) is analytic 
at {Fx{s()), Sq), hence we can apply Fact 5.2 again to conclude that Fx is analytic at sq as 
well. □ 

The next lemma ensures that we can ignore certain degenerate cases appearing in the 
course of the analysis of (3.5). It is the analytic interpretation of the probabilistic results in 
Section 2. 

Lemma 6.3. 1 is a singular point of a{s). If c = cq, then a"{l) = +00. 

Proof. If c = Co, the second assertion follows from Theorem 1.1 or from Neveu's result that 
E[Zx log^ Zx] = 00 for X > (see the remark before Theorem 1.1). This implies that 1 
is a singular point of a{s). If c > cq. Proposition 2.1 implies that E[s^^] = 00 for every 
s > 1, whence 1 is a singular point of the generating function Fx{s) by Pringsheim's theorem 
([13], Theorem IV. 6, p. 240). By Corollary 5.9, it follows that 1 is a singular point of a(s) as 
well. □ 

The next theorem is the core of the proof of Theorem 1.2. 
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Theorem 6.4. Under the assumptions of Theorem 1.2, for every 93 G (0, vr) there exists 
r > 0, such that a{s) possesses an analytical extension (denoted by a{s) as well) to G{ip,r). 
Moreover, as 1 — s — t- 1 in G{ip,r), the following holds. 



If d = 1, then 



_s log log ^ 

■logi ''°'(logi)2 ' - \{lo^ 



a{l -s) = -cos + co^ - cos^^-f-^ + O —-j- . (6.2) 



If d > 1, then there is a K = K{c,f) G C\{0} and a polynomial P{s) = Yl\i=2'^'ns"' , 
such that 

i/d^N: a{l-s) = -\cS + P{s) + Ks'^ + o{s'^), (6.3) 

ifdeN: a{l- s) = -XcS + P{s) + Ks'^logs + o{s'^). (6.4) 



Proof of Theorem 6.4- We set b{s) = a(l — s). By (6.1), 

-6'(s)6(s) = 2c6(s) + 2(1-5-/(1 -s)) onD(l,l). (6.5) 

Since / is analytic at 1 by hypothesis, there exists < ei < 1 — q and a function g analytic 
on D(0,ei) with g{0) = g'{0) = 0, such that /(I - s) = 1 - (m + l)s + g{s) for s G D(0,ei). 

As a first step, we analyse (6.5) for real non-negative s. Since ei < 1 — g', 6(s) < on 
(0,ei), whence we can divide both sides by b{s) to obtain 

db -2cb - els + 2g{s) 

— = ^ on(0,ei). (6.6) 

Introduce the parameter t{s) = f^^ z^pji ^ ^ (0,£i], such that t(ei) = 0, t(0+) = +00 
and t(s) is strictly decreasing on (0, ei]. There exists then an inverse s{t) on [0, 00), which 
satisfies s'{t) = b{s{t)). Hence, we have 

db db ds „ w X 9 / N / / NX 

di " d^ dt " ~ " + ^^'^^^^ ^ ' 

In matrix form, this becomes 

iC)-<:)-(T). -^(-f-o^). (-) 

for t G (0,00). Note that this extends (4.2) to the subcritical case. This time, the Jordan 
decomposition of M is given by 

A-H,A=(f _\), A^(-f -^A. .foe. (6.8) 



and by (4.3), if c = cq. Setting 
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transforms (6.7) into 



^ = -XcB + [B,S]2, ^ = -X,S+[B,S]2, ifoco, (6.10) 
^ = -coB + S+[B,S]2, ^ = -coS+[B,S]2, if c = cq, (6.11) 
for t € (0,cx3). Furthermore, by (6.9), we have 

s = B + S, (6.12) 

^^UX-XcrHb + Xs), ifd>l, ^g^^^ 
1 6 + cqs, if d = 1, 

^^UXc-XyHb + X,s), ifd>l, 
\-b+{l-co)s, if d = 1. 

From now on, let £2, £3, . • • be positive numbers that are as small as necessary. By the strict 
convexity of b and the fact that b'(0) = —Ac by Lemma 2.4, equation (6.13) implies that S is 
a strictly convex non-negative function of s on [0,£2). This implies that the inverse s = s{S) 
exists and is non-negative and strictly concave on [0, £3). It follows that t{S) = t{s{S)) exists 
on [0, £4). Equations (6.10) and (6.11) then yield for S £ (0,£4), 

dB dB + [B,S]2 ^ 

dS= S+[B,S]2 ' ^'-''^ 
dB B-c^'S+[B,S]2 

dS= S+[B,S]2 ' = ^'-''^ 

By (6.12) and the fact that s{S) is strictly concave, S is a strictly concave function of S 
as well, hence strictly monotone on (0, £5). We claim that B{S)'^ = o{S) as S" — )■ 0. For d > 1, 
one checks by (6.13) that S{s) ~ s, as s — > 0, whence B{S) = o{S), as S* — 0, by (6.12). 
If d = 1, then b'{0) = — cq by Lemma 2.4 and 6"(0) = +00 by Lemma 6.3. Equation (6.13) 
then implies that S{s)/s'^ — )• +00 as s — )• 0, whence s{S) = o(^/S). The claim now follows by 
(6.12). 

Proposition A. 4 now tells us that there exists a function h{z) = [z]2, such that the function 
s{S) = S — h{B{S)) has an inverse ^(s) on (0,£6) and 6(5) = i?(5(s)) satisfies the Briot- 
Bouquet equation 

^^,^U + [M]2, ifd>i, 

- Cq 5+ [b,5]2, if d = 1, 

on (0,£6). By Fact 5.4 and Remark 5.5, there exists then a function tp{z,u) = u + rz + [z,'ul\2, 
r G C, such that b(s) = '(/'(s, u(s)), where 

u{z) = Cz'^, if d ^ N and u{z) = Cz'^ log z, if d G N, 

for some constant C = C{c, f) £ C (the form of u in the case d G N can be obtained from the 
one in Fact 5.4 by changing ^, C and K). Moreover, comparing the coefficient of 5 on both 
sides of (6.17), we get, if d > 1, r = dr, whence r = and ifd=l: r + C = r — Cq ^, whence 
C = -c^\ 
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Assume now d > I. Then b = u{5) + [s,ti(s)]2. Recall that B = b and S = 5 + h{b). By 
(6.12), 

s = B + S = b+5 + h{b) =5 + n(s) + [s,ti(s)]2, 

such that s'(s) = 1 + o(l) and s(s) = s + [s]2 + 0(5'^), as s — > 0, where 7 = (d + [dJ)/2, if 
d and j = d — 1/2, if d G N. By Lemmas A. 6 and A. 8, for every € (0, vr) there exists 
ro > 0, such that the inverse s(s) exists and is analytic on H{ipQ,rQ) and satisfies 

s(s) = s + [s]2 + o(s'^), as s — ;> 0. 

This entails that 

n(s) = Cs'^ = C{s + o{s)f = Cs'^ + o{s'^), if d ^ N, 

n(s) = C5'^\og5 = C{s + o{s^l'^)f\og{s + o{s)) = C s'^ log s + o{s'^) , if dGN\{l}, 

s" = [s\2 + o(s^+i) + o(s^') = [s\2 + o{s'^), for all n>2. 

It follows that 

b(s) = b(s(s)) = 'u(s) + [s]2 + o{s'^), as s 0. 
We finally get by (6.9), 

b = -XB - X,S = -Xcs + {X,-X)b = ~XcS + (Ae - K)u{s) + [s]2 + o{s'^), 

which proves (6.3) and (6.4). 

If d = 1, recall that u{z) = CQ^^zXog^ and b = n(s) + rs + [s,u(s)]2 for some r G C. By 
(6.12), 

s = B + S=b+5 + h{b) = u{5) + (r + l)s + [s, ti(s)]2, 

such that s'(s) = Cq^ log(j) + 0(1) and s(s) = Cq ^slog j + (r + 1)5 + o(s). Lemma A. 6 now 
implies that for every (pQ € (0, vr) there exists ro > 0, such that the inverse 5{s) exists and is 
analytic on H{ipQ,rQ). Now, by (6.9), 

b = -Cos + S = -Cos + s + h{b) = -C0S + 5 + 0(s^/^). 

Lemma A. 9 now yields (6.2). □ 

Remark 6.5. The reason why we cannot explicitly determine the constant K in Theorem 6.4 
is that we are analysing (3.5) only locally around the point 1. Since the solution of (3.5) 
with boundary conditions a{q) = a(l) = is unique (this follows from the uniqueness of the 
travelling wave solutions to the FKPP equation), a global analysis of this equation should be 
able to exhibit the value of K. But it is probably easier to refine the probabilistic arguments 
of Section 2, which already give a lower bound that can be easily made explicit. 

The asymptotics established in Theorem 6.4 for the infinitesimal generating function can 
now be readily transferred to the generating functions Fx{s). 

Corollary 6.6. Under the assumptions of Theorem 1.2, for every x > and ip G (0,7r) there 
exists r > 0, such that Fx{s) = E[s'^'^] can be analytically extended to G{ip,r). Furthermore, 
the following holds as 1 — s — t- 1 in G{ip,r). 
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If d = 1, then 

s log log ^ 

logi (logi)2 / ■ - \ (logi) 



-s) = l- e^'^^s + coxe^^^ ( ^ - +0 — ^ . (6.18) 



— If d > 1, then there is a polynomial Px{s) = X]i,=2Cn'S", such that 

ifd(^n: - s) = 1 -e^=^s + P,(s) + K^(is'^ + o(s'^), (6.19) 

ifd^N: Fx{l - s) = 1 -e^-''s + Px{s) + Kxs'^ log s + ois'^), (6.20) 

where = K{e^''^ — e'^''^)/(Ac — Ac), with K being the constant from Theorem 6.4- 

Proof. Let < ifo < if. By Theorem 6.4, there exists tq > 0, such that a{s) can be analytically 
extended to G((^Oi^o) ^ind satisfies the hypothesis of Corollary 5.10. It follows that there 
exists r > 0, such that Fx{s) can be analytically extended to G{if,r) and maps G{if,r) into 
G{ipo,rQ). Hence, the functions 

f*{l-r)dr, 

where f*{s) is defined as in (5.4), are analytic in H{'jr — (p,r). In what follows, we always 
assume that s G H{7r — ip,r). Appearance of the symbols ~,0,0,o means that we let s go 
to in H(7r — ip, r). 

First of all, we note that by Proposition 5.8, we have 

w{s) = se^=^ exp(/(s)) = se^-^ (l + I{s) + ^ . (6.21) 

k=2 

Now assume d> \. By Theorem 6.4, a(l — s) = —\cS + [s]2 + u{s) + o{s'^), where u{s) = Ks'^ 
or u{s) = Ks'^logs, according to whether (i ^ N or d G N, respectively. It follows that 

Now, f^^^^ o{r'^~'^) dr = o{s^~^), since w{s) ~ se^"^ by Lemma 2.4. Thus, 

I{s)= r(l-r)dr=Ks),s]i- / ^dr + o(s'^-i). (6.22) 

Js Js ^cf 

Since j"^^^^ r~'^u{r) dr = 0(s'^~^ log s), equations (6.21) and (6.22) now give 

w{s) = se^=" ^1 + [w{s), s]i - J^^'^ + o{s''-^)^ . (6.23) 

If d > 2, we deduce that w{s) = se^"^ + o(s'^/^). Straightforward calculus now shows that 

/ >^ dr = + [w{s\ s], + o{s'~% (6.24) 
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and (6.23) and (6.24) now yield 

w{s) = se^'^'^ + [w{s), s]2 - K^u{s) + o{s'^). 

Repeated application of this equation shows that w{s) = se^"^ + [s]2 — Kxu{s) + o{s'^), which 
yields (6.19) and (6.20). 

In the critical case d = 1, Theorem 6.4 tells us that 



s \ logi (logi)2 \ (log-i^2 ; ; 



Write A = Ac = cq. For our first approximation of w{s), we note that 



sc^- 1 , 1 r^^^ 1 , Xx 



I{s) — Y dr ~ Y / - dr — , , 

Js r log i log ^ Js r log i 

hence, by (6.21), 

w{s) = se'^ (l-^+o f I I . (6.26) 



log- \log- ^ 
To obtain a finer approximation, we decompose I{s) into 

-r)dr+ r (1 - r) dr =: h{s) + his). 

We then have 

log 7 (logf)^ y(iogi)^y 

and, because of (6.26), 

-h[s) ~ / — r dr 



rlogi (logi)2' 

Plugging this back into (6.21) finishes the proof. □ 

Proof of Theorem 1.2. Let x > 0. We want to apply the methods from singularity analysis 
reviewed in Section 5.3 to the functions a and Fx, if (5 = 1, or the functions h and from 
Lemma 6.2, \i 8 > 2. Let ip G (0, 7r/2). By Theorem 6.4 and Corollary 6.6, there exists 
ro > 0, such that a and can be analytically extended to G{ip,rQ), which implies that 
for some ipi € (93, 7r/2) and ri G (0,r), /i and can be extended to G{ipi,ri), as well. 
Moreover, by Lemma 6.2, each of these functions is analytic in a neighbourhood of every 
point of C = {s € dD : |1 — s| > which is a compact set. Hence, there exists a finite 

number of neighbourhoods which cover C. It is then easy to show that there exists r > 0, 
such that the functions are analytic in A{ipi,r). 

If (5 = 1, we can then immediately apply Facts 5.6 and 5.7, together with the asymptotics 
on a and F^ established in Theorem 6.4 and Corollary 6.6, to prove Theorem 1.2. 

If 5 > 2, let q{s) be the inverse of s 1— >• s"^ in a neighbourhood of 1, then h{s) = a{q{s)) / q{s) 
near 1, by Lemma 6.2. But since q'{l) = 1/5, we have 

h{l -s) = a{l - ds + C2S^ + C3S^ + • • • ))(1 + c'lS + c'^s + ■ ■ ■ ), 
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for some constants c,n,c^, and so equations (6.2), (6.3) and (6.4) transfer to h with the 
coefficient of the main singular term divided by 6'^. We can therefore use Facts 5.6 and 5.7 
for the function h to obtain the asymptotic for {psn+i)neN in Theorem 1.2. In the same way, 
equations (6.18), (6.19) and (6.20) yield asymptotics for hx, such that we can use again Facts 
5.6 and 5.7 to prove the second part of Theorem 1.2. □ 



A Appendix 



A.l A renewal argument for branching diffusions 

Let W 



{Wt)t>o be a diffusion on an interval with endpoints a' < < a, such that 
lima;4olP''[7o < t] = 1 for every t > 0, where Tq = inf{t > : Wt = 0} and Wq = x, 
P^-almost surely. For x € (0, a), and only in the scope of this section, we define to be the 
law of the branching diffusion starting with a single particle at position x where the particles 
move according to the diffusion W and branch with rate /3 according to the reproduction law 
with generating function f{s). Moreover, particles hitting the point are absorbed at that 
point. Denote by Z the number of particles absorbed during the lifetime of the process and 
define Usix) = P^ [s^] for s £ [0, 1) and x £ (0, a). 

Lemma A.l. Let s € [0, 1) and Q he the generator of the diffusion W . Then 

Qus = f3{us — f o Us) on (0, a), with Us{0+) = s. 

Proof. The proof proceeds by a renewal argument similar to the one in [22]. As for the BBM, 
for an individual u, we denote by its time of death, Xu{t) its position at time t and the 
number of u's children. Define the event A = {3t G [0, C0) '■ ^zit) = 0}. For s G [0, 1) we 
have by the strong branching property 



Usix) 



sP'^{A)+E'' 



where W = {Wt)t>o is a diffusion with generator Q starting at x under P^', Tq = inf{t > : 
Wt = 0} and the random variable ^ is exponentially distributed with rate (3 and independent 
from W. Setting v{x) = ¥^{Tq < ^) we get by integration by parts 



v{x) 



/3e-^*P^(ro < t) dt 



■(To G dt) = 



-/3To 



and therefore Qv = (3v on (0,a) ([8], Paragraph II. 1.10, p. 18). 

Denote the /3-resolvent of the diffusion by Rj^. By the strong Markov property. 



mfius{Wii)),^<To]=E^ 



^e-^'f{us{Wt))dt 



To 



l3e-f''f{usiWt))dt 

= (3Rp{f o us){x) - l3E^\e-^^'^]Rp{f o ^x,)(0), 
hence Us = Cs^pv + [iRp{f o Us), with Cs^(3 = s — fiRj3{f o Us){0). It follows that 

GUs = f3Cs,i3V + P'^Rpif ° Us) - I3{f o Us) = P{us - fous) ou (0,a). 
By the above hypothesis on W, P^(To < .^) — )■ 1 as x ^ 0, whence Us{0+) = s. 



□ 
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A.2 Addendum to the proof of Theorem 1.1 

With the notation used in the proof of Theorem 1.1, recah that for some constant K > we 
have 

{(p' + co(l)){x) ~ (p{x)/x ~ Ke~^°^, as X ^ oo. 

In what follows, formulae containing the symbols ~ and o() are meant to hold as s J, 0. The 
above equation yields 

a(l -s) = = -cos + + com-\s)) = -cos + ^"^^ + ""^^^ . (A.l) 

Now, by (3.5), we have 

a'(l — s)a{l — s) = 2coa(l — s) + CqS — g{s), 

where we recall that g{s) was defined as g{s) = 2(/(l — s) — 1 + /'(l)s). From the above 
equation, one gets 

a"{l -s) = -(a(l - s))-^[{coa{l - s) + - g{s))' - g'(s)a(l - s)^), 

and an application of Lemma 4.1 and (A.l) yields (4.7). 

Kolmogorov's forward and backward equations (3.2) and (3.3) give 

a{s) ' 

and taking the derivative on both sides of this equation gives 

F^{s) = ^^^ia'{FAs))-a'{s)). (A.2) 
By (4.7) and F^(l) = E[Z^] = e'^"^, we get 



-1-Fa;(l-S) 

a\F^{l - s)) - a'{l -s) = - I a"(l - r) dr ~ 



Cr,X 



(log S)2 ■ 

This equation, together with (A.2) now yields 



F^il - s) 



-CQe'^^^s I CnX 



CqS^ y (logs)2y ' 
which is (4.8). 

A. 3 Reduction to Briot Bouquet equations 

In this section, we show how one can reduce differential equations as those obtained in the 
proof of Theorem 1.2 to the canonical form (5.1). It is mostly based on pp. 64 and 65 of [5]. 

Lemma A.2. Let A G (0, 1] and p € C. Then the equation 

w' = + . (A.3) 

z + pw + [w, Z\2 

has an analytic solution w{z) = [z]2 in a neighbourhood of the origin. 
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Proof. We choose the ansatz w = z ■ wi. This transforms (A. 3) into 

, Xzwi + z'^[wi,z]o Xwi + z[wi,z]o 

ZWi +Wl = JTT 7- = ; f . 

z + pzwi + Z^[Wl, z\o l + [wi,z\i 
Writing the inverse of the denominator as a power series in wi and z, this equals 

{Xwi + z[wi,z]o){l + [wi,z]i) = Xwi +rz + [wi,z]2, 

for some r G C. This finally yields 

zw'i = (A - l)wi + rz + [wi,z]2. 

Since A — 1 is not a positive integer, this equation now has an analytic solution wi{z) = [z]i 
by Fact 5.3, whence w{z) = zwi{z) = [z]2 solves (A. 3). □ 

Remark A. 3. The important point in Lemma A. 2 is that the coefficient of z in the numerator 
of (A. 3) is 0, which is why w'{z) = 0. 

Proposition A. 4. Let A > 1 and p £ M. Suppose w{z) is a strictly monotone real-valued 
function on (0, e), e > 0, with w{z)'^ = o{z) as z — )• and satisfying 

I Xw+pz^\w,z\2 , . , . 

w = -. ^ on (0,e). (A. 4) 

Then there exists h{z) = [z]2 and Si > 0, such that I = z — h{w) has an inverse z = z{^) on 
(0, ei) and such that 

3— = Aw; + p3 + [w,3]2 on(0,ei). (A.5) 

da 

Proof. By hypothesis, w{z) is monotone on (0,e) and therefore possesses an inverse z = z{w) 
on (0,5), 5 > 0, which satisfies 

dz X-^z + [w,z]2 , . . 

~r~ = ^~^ F T- (A-d) 

aw w + pX '-z + [w,z\2 

By Lemma A. 2, there exists then an analytic solution z = g{w) = [w\2 to (A. 6), since 
A~^ G (0, 1] by hypothesis. Setting i = z — g{w) transforms (A. 6) into a differential equation, 
which has 3 = as a solution, hence it is of the form 

d3 A-^3 + 3[w,3]i 



dw w + pX '^i + [w,])\2 ' 

We have d3/dz = 1 + g' {w{z))w' {z) = 1 + 0{w{z)w' {z)). By (A.4), 

w{z)w' {z) = 0{w{z)'^ /z + w{z)) = 0(1), 

by hypothesis. Hence, there exists ei > 0, such that i{z) is strictly increasing on (0, ei) and 
therefore has an inverse. Thus, w{'^) = w{z{i)) satisfies 

du; w + pX~'^i + [w,i\2 . 

-5- = — rrrrm on(0,e2), 

d3 A 13(1 + [w,3]i) 

for some £2 > 0. Expanding (1 + [tL',3]i)~^ as a power series at (ii),3) = (0,0) gives (A.5). □ 
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A. 4 Inversion of some analytic functions 

The results in this section are needed in the proofs of Corollary 5.10 and Theorem 1.2. 

Lemma A. 5. Let 99 G (0, vr), r > and h be an analytic function on H{ip, r) with h{z) = o{z) 
as z — )• 0. Then there exists ri > 0, such that for all zi,Z2 S Hif, ri), 

log zi - log 22 + / h{z)dzy^O. 

J Z2 

Proof. Let zi,Z2 G H{{p,r). Write Zi = with > 0, (pi £ (—93, 99), i = 1,2. Define the 

paths 

= a2e'(*^i+(i-*)'^2) and 72(i) = (toi + (1 - t)a2)e''^i , tG[0,l], 
such that their concatenation forms a path from Z2 to zi in H{ip,r). Then 

/ /i(s) ds = — 932I • ^20(1/02) and / /i(s) ds = | log ai — log a2|o(l). 

As a consequence, 



h{s) ds 



Z2 



- <^2 



+ I logai - loga2|)o(l) < \/2| logzi - log2;2|o(l). 



This proves the statement. □ 

Lemma A. 6. Let r > and ip € (0, vr]. Let g and h he analytic functions on H{{p,r) with 
g'{z) = 1 + 0(1), h\z) = logi + 0(1), g{z) and h{z) ^ as z ^ in H{ip,r). Then 
for each cpo G (0, and ipi £ (ifo, ip) there exist ro, ri > 0, such that g and h are infective on 
H{(pi,ri) and the images of H{ipi,ri) by g and h contain H{ipQ,rQ). 

Proof. By hypothesis, g{z) = z + o{z) as 2; — )• in H{ip,r), whence a,rgg{z) = argz + o(l). 
Thus, there exists ri > 0, such that g{H{(pi,ri)) C C\(— oo,0]. 

Suppose that there exist zi, Z2 G H{(pi,ri), such that g{zi) = g{z2). Let 7 be a path from 
Z2 to zi in H{ipi,ri). Then (7 o 7 is a loop in C\(— cxd,0], whence 

= [ -dz= [ ^dz = logzi-logZ2+ / o(i)dz. 

By Lemma A. 5, we can choose ri so small, that this equality cannot hold, whence g is injective 
on H{ipi,ri). 

Since g{z) — )• and avg g{z) = argz + o(l) as 2; — )• 0, there exists tq > 0, such that 
g{dH{ipi,ri)) encloses H{ipo,rQ). Now, since g is injective on H{ipi,ri), H{ipi,ri) and 
g{H{ipi,ri)) are conformally equivalent, whence g{H{ipi,ri)) is simply connected. It follows 
that g{H{ipi,ri)) D H{ipo,ro). 

Exactly the same arguments hold for h, since h{z) = z(log -+0(1)) by hypothesis, whence 
avg h{z) = argz + o(l) and h'{z)/h{z) = l/z + o(l/z) as z — )• 0. □ 

Lemma A. 7. Let r > 0, (/? € (0,vr] and t G M. Let g be an analytic function on H{ip,r) with 

g'{z) = - + —^ ( ) , as z in H{(p,r), 

z 2: log z \z\ log zp J 

for some c G M and 7 > 1. Then for each < pQ < ipi < ip there exist rg, ri > 0, such that g 
is injective on H{ipi,ri) and g[z) + t G g{H{ipi,ri)) for every z G H{ipQ,rQ). 
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Proof. By the hypothesis on g, we have for zi,Z2 G H{ip,r), 

g{zi) - g{z2) = log zi - log Z2 + o{l/z) dz. 

By Lemma A. 5, there exists therefore ri > 0, such that g is injective on H{ipi,ri). 
Since l/(x| logx|'^) is integrable near 0, we have 

g{z) = logz + clog(log i) + o(l), as z 0, 

where we assume without loss of generalisation that the constant of integration is 0. It follows 
that Jieg{z) — >■ — oo and Im.g(z) = argz + o(l) as z — ?• 0, since c € M. Hence, there exists an 
i? S M, such that g{dH{ipi,ri)) encloses the strip 5 = S-{R,ipi). As in the proof of Lemma 
A. 6, it follows that S C g{H{ipi,ri)). Furthermore, again by the asymptotics of Re^ and 
Im(jr, there exists tq > 0, such that g{s) +t G S for every s € G((/9o>'"o)- This concludes the 
proof. □ 

Lemma A. 8. Let w{z) be an analytic function on an open subset o/C\(— oo,0], such that 
w{z) — )• as 2: — )• and 

z = w + a2uP' + • • • + Onw"^ + o(u;^), as z — )• 0, 
for some n e N, 7 > n and 02, . . . , G C. Then there exist 62, . . . , 6„ G C, such that 

w{z) = z + b2z'^ -\ l-buz"" + o{z'^), asz^O. 

Proof. For every i S N, we have by hypothesis 

z' = w' + ai^i+iw'^'^ H h ai^nw"" + o{w'^), 

for some flj j^i, . . . , ai,n £ C. For 2 ^ /c ^ 71, define recursively (with bi — 1) 

bk = -{ai,k + b2a2,k H h 6fc-iafc-i,A:)- 

Then, z + 62-2^ + • • • + b"'z^ = w + o(w;'^). The statement now follows from the fact that 
w{z) ~ 2; as z — )• by hypothesis, whence o{w'^) = o{z"'). □ 

Lemma A. 9. Let w{z) be an analytic function on an open subset o/C\(— oo,0], such that 
w{z) — )• as z — )• and 

cz = w log ^ + Cw + o{w) , as z — )• 0, 
for some constants c > 0, C G C. Then 



cz 



, iogiogi + (:;-iogc + o(i) \ 

w = Y 1 1 ) as z — ^ 0. 

log I V logi / 



Proof. Set f{z) = w{z)/z. By hypothesis, logz ^ logu; = logz + log/(z), whence log/(z) 
o(log2;). Now define g{z) by 

wiz) = Y9iz), 

log 7 
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such that log g{z) = log f{z) — log log ^ = o{logz). By hypothesis, 

cz ~ wlog i = 7^^s{z) (log log i + log i - logc - log 5^(2;)) ~ czg{z), 



cz = -^g{z) ( log log i - log cz + C + o(l) ) , 



whence g{z) ~ 1, which implies log g{z) = o(l). It now follows from the hypothesis that 

cz 
log-^ 

whence 

/ loglogi + C-logc + o(l)\"^ 

^'^*=r — — j ^ 

The statement now follows from the series representation of (1 + z)~^ at z = 0. □ 
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